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ABSTRACT 


Blanket chemical vapour deposited tungsten (W-CVD) is routinely employed to 
fabricate reliable contacts for sub-micron multilevel metallization. The deposition is 
achieved by hydrogen reduction of tungsten hexafluoride. In case of a trench, as an 
example, as the aspect ratio, defined as the ratio of depth to width of a trench, increases, a 
void may form, leading to serious problems in the reliability of integrated circuits (ICs). 

In the present work we investigate a novel idea of imposing a controlled thermal 
gradient along the vertical direction of the trench, in a manner that the base of the trench 
is hotter than its mouth. 

In the mass transport and chemical reaction model presented here, Knudsen flow, • 
viscous flow and ordinary diffusion are all taken into account. The mass flux in the two- 
dimensional model is expressed according to the “dusty gas” model. The model 
equations are numerically solved via a finite difference discretized approach for the 
parameters corresponding to tungsten deposition by hydrogen reduction of tungsten 
hexafluoride. It also accoimts for the changes in the trench shape that occur during the 
deposition of tungsten on the trench boundaries. 

The isothermal deposition profile computed for trench of different aspect ratios 
show the formation of a void at deposition temperatures 673 and 773 K. In order to attain 
void firee deposition, a linear temperature gradient is imposed keeping the temperature at 
mouth of the trenches at 673 and 773 K, respectively. At the most the base of the trench 
is 3 K hotter. We find excellent step coverage for different aspect ratios in the range 2-5 
(ratio of trench depth to width of the trench), at deposition temperatures of 673 K. Even 



after an application of thermal gradient, void still forms at 773 K, while at 573 K, also 
investigated in this work, no void forms even under isothennal process. 

We also discuss the dielectric material that could support the thermal gradient 
applied in this work. We suggest that either commonly used Si02 layer will have to be 
processed differently, or new material will be used to account for required value of 
thermal conductivity. In this regard, polyamide or fluoropolyamide dielectrics may be 
more suited, which are also candidates for interlayer dielectric for their dielectric 
properties. 

Keywords: CVD; contacts; step coverage; ICs 
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CHAPTER 1 


INTRQDIJCTTON 


1.1 BACKGROUND 

The formation and behavior of metal-semiconductor contacts play a vital role in the 
very large scale integrated (VLSI) technology. For metal oxide semiconductor (MOS) 
transistors, metals and their silicides are used as gate electrodes, contacts, diffusion 
barriers and in intercoimect metallization. 

With advances in integrated circuit technology, device dimensions are being 
scaled down to submicron levels. Also the integrated circuits (ICs) require chips having 
fast switching speed, high reliability and low cost per function. The larger size and high 
complexity chips require closely placed long interconnect lines with a small area of cross 
section. As a result, the resistance-capacitance (RC) time delay becomes significant. The 
RC time constant of a interconnect is given by: 

RC = RsL^s/t 

where Rs and L are the sheet resistance and length of the connection, respectively; t and e 
are the thickness and dielectric constant of the insulating film, respectively. Thus the 
performance of a large circuit will be seriously affected by the interconnections present in 
it. In order to achieve higher switching speed the sheet resistance (Rs) of the interconnect 
should be reduced. 
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Heavily doped polycrystalline silicon (pKily-Si), aluminum (Al), and aluminum 
alloys have been widely used in VLSI circuits. However, the high electrical resistivity 
value (p = 500 |a,ohm.cm) of poly-Si results in considerable signal propagation delay and 
joule heating. This limits its suitability for metallization of submicron size devices. In the 
case of Al, in- spite of its low electrical resistivity (p = 2.7 pohm-cm), it also has two 
major drawbacks which limits its use in ICs. First, at high current density, 
electromigration destroys the physical integrity of aluminum interconnect lines. Second, 
Si has a tendency to diffuse into aluminum grain boundaries, resulting in spikes of 
aluminiun. These Al spikes cause leakage and possible shorting of the junctions[l]. 

Another severe limitation of using Al is its nonconformal coverage over the 
interconnects. Because of the inherent conformal coverage properties of chemical vapour 
deposition (CVD), hmgsten deposited by CVD is a candidate process in interconnects 
having high aspect ratios. The following features of tungsten (W) make it attractive for 
VLSI circuits. 

• Deposition can be selective or non-selective. 

• It is easy to pattern non-selective tungsten films by wet or dry etching techniques. 

• Excellent conformal coverage is shown by non-selective or blanket depositions 
process. 

• Films can be deposited in hot wall or cold wall CVD reactors, over a wide range of 
temperatures. 

• Excellent electromigration resistance of W films compared to Al. 

• Good compromise between prevention of jimction spiking and an increased electrical 
resistivity (compared to Al). 
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• The thermal expansion coefficients of tungsten and silicon are similar, which reduces 
the stress at the interface. 

Tungsten hexafluoride (WFe ) is a suitable precursor gas for tungsten 
deposition in ICs. CVD of tungsten can be achieved both by silicon reduction and 
hydrogen reduction of WFe, according to the following reactions: 

2WF6 (g) + 3Si (s) = 2W (s) + 3SiF4 (g) 

WFe (g) + 3H2 (g) = W (s) + 6HF (g) 

According to the free energy data, silicon reduction of WFg is thermodynamically 

favourable over the hydrogen reduction reaction. However, the silicon reduction reaction 

0 

is self-limiting, that is, after a film of about 100-200 A , the reaction stops. Subsequently, 
thicker tungsten films can be grovra on silicon substrates by Ha reduction of WFe. There 
is no inherent limit on the thickness of tungsten films that can be deposited by H 2 
reduction reaction. 

Tungsten deposition in the interconnects of ICs can be achieved selectively or 
non-selectively. The non-selective deposition is often termed as blanket deposition. In 
blanket tungsten deposition, after metal deposition of W by silicon reduction reaction, 
tungsten is deposited non-selectively on the already deposited metal in the trench as well 
as on the dielectric by hydrogen reduction reaction. The excess tungsten deposited on the 
plug and the other areas is etched off thereafter. On the other hand, in case of selective 
deposition, tungsten is deposited only in the trench and not on the dielectric or other, 
insulators. This selective deposition could be used to eliminate mask and etch steps, 
required in the case of blanket deposition, to remove excess metal. Here tungsten 
deposition first proceeds by silicon reduction and subsequent W deposition occurs on the 
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already grown tungsten by hydrogen reduction. However, this process requires specific 
processing conditions and a tight control of process variables. Moreover, the selective 
nature of tungsten deposition often breaks down before the deposition is complete. 

Thus frequent loss of selectivity makes the selective tungsten deposition process 
less popular in commercial reactors and, the blanket tungsten deposition is the preferred 
approach. 

1.2 LITERATURE REVIEW 

In the following, a review of literature is provided on tungsten CVD for filling 
intercoimects in VLSI circuits, concentrating especially on a void that is formed inside 
the trench while depositing tungsten from its precursor gas WFe. 

In an integrated tungsten deposition and etchback plug fabrication process, in a 
cold wall single wafer CVD chamber[2], at a temperature above 300°C, it is possible to 
obtain uniform tungsten plugs in devices without degrading the electrical performances 
of the contacts (compared to standard aluminum case). 

Tungsten deposition is usually carried out by CVD process, which in comparison 
to PVD, yields a better step coverage. The deposition of tungsten starts by reduction of 
WFe with silicon substrate. But due to the self-limiting nature of this reaction, 
subsequent tungsten deposition is carried out by H 2 reduction. 

Tungsten deposition can be done in two ways: blanket and selective. A review of 
the current W plug process is provided by Ireland[3], who discusses the role of dopants in 
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obtaining good metal/silicon interface, the associated cleans, barrier layers and then 
tungsten deposition and etch back processes. 

The selective deposition has received much attention in the past. But due to the 
tight process conditions[4], reactor design and selectivity, this process has not been fully 
commercialized yet. On the other hand, blanket deposition which suffers from the 
drawbacks of incomplete filling, is still practised in majority of VLSI fabrication 
processes. 

Pauleau and Lami[5] investigated the kinetics and mechanism of W deposition 
both by silicon and hydrogen reduction. Under silicon reduction it was found that[5,6], 
the deposition is purely selective and no encroachment of W film occurs along Si/SiOa 

interface. It was also found that within a temperature range of 300-475°C; the film 

0 

thickness is limited between 100 and 150 A, indicating self-limiting nature, of silicon 
reduction. Another related work is proceeded by Yeh et al.[7], in which they were able to 
deposit selective W film by silane reduction of tungsten hexafluoride. The selectivity 
becomes excellent at a deposition temperature of 280-350°C and a flow rate of SiEUAVFg, 
less than 0.6. However, at temperatures above 350°C, W starts nucleating on Si02 
surface, because of the by-product SiF 2 with WFe, leading to selectivity loss. In WF 6 -H 2 
deposition, the rate limiting mechanism was foimd to be dissociation of hydrogen on 
substrate and the reaction is one-half order in H 2 concentration and zero order in WFa 
concentration[8]. But this rate expression does not accoimt for the fact that rate must go 
to zero as the concentration of WF6 goes to zero. Therefore, Jain et aL[9] modified the 
rate expression as; 
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Rate = 0.00827 exp 


f 8800^ 

Ph' 


1 T ; 


PwFg (1 + kfPwFfi 


mol/cm /s 


( 1 ) 


where T is the temperature in FI, Ph and are the hydrogen and tungsten hexafluoride 


partial pressures in Torr, respectively, P^p^is the partial pressure of tungsten hexafluoride 
at the top of the wafer and kf is a constant whose value is 1 000 Torr'* . 

A low temperature (250-500°C) is required for selective deposition of W by H 2 
reduction. Moreover, the concentration of both H 2 and WF6 must be low, that is, 
deposition rate must be slow. But in this case the selectivity is maintained only for a 
short deposition time (15-30min). Many investigations tried to find out the reason for 
loss of selectivity. Pauleau and Lami[5] associate loss of selectivity to presence of high 
concentration of by-product gas HF, whose interaction with SiOz results in nucleation of 
W on oxide. Broadbent and Stacy[6] finds that the loss of selectivity is mainly due to 
presence of silicon containing by products, such as SiF 4 . The reason for loss of 
selectivity is also explained by Kuijlaars et al.[10] and Werner et al.[ll]. Experimental 
evidence reveals that nucleation of tungsten occurs in the areas farther from the trench, 
which means the loss of selectivity is primarily due to gas phase diffusion rather than 
surface diffusion. The nucleation accelerates by the presence of additional W surfaces in 
the reactor[ll]. 

Here it will be useful to cite two theories that explain the loss of selectivity. In 
the first, HF is thought to be responsible, which reacts with Si02 to form SiFx (x<3) over 
a tungsten catalyst. This is followed by the interaction of reactive SiFx with WFc, to form 
a deficient WFy, which when adsorbed on oxide surface form tungsten islands. However, 
experiments do not match well with the degree of selectivity loss. Another theory, which 
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deals with observation of adsorbed species on oxide, suggests presence of sub-fluorides 
WF5[12]. 

Apart from the experimental studies, some modeling of hmgsten CVD is also 
reviewed. Ulacia et al.[13] studied the physical and chemical simulation of a tungsten 
CVD reactor using fluid dynamics equations. The effect of thermal diffusion has been 
completely neglected. This study provides a better understanding of the CVD process, 
optimizing the variables that affect deposition rate and film uniformity. A similar work is 
carried out by Cale et al.[14], where they have decreased the processing time for a given 
step coverage by var 3 dng the deposition rate in a prescribed manner during the process. 
A theoretical model by Joshi et al.[15] illustrates that the growth rate depends on the 
trench size, aspect ratio and WFe concentration. 

In another work, Kuijlaars et al.[16] developed a model to study the influence of 
process condition and the concentration distribution of tungsten sub-fluorides, which are 
believed to be the reasons for loss of selectivity. They found that low temperature, low 
pressure and large total gas flow rate enhances selective deposition. However, their study 
left some unresolved problem regarding the chemistry of the species, to which some ad 
hoc modifications were made. 

From several modeling and experimental studies on selective W deposition, it is 
obvious that tight process control is needed for enhancing selectivity e.g. 

• Low pressure followed by low growth rates. 

• Low temperature. 

• Cooled reactor wall to minimize reaction and hence product gases. 

• Coating the back of the wafer for avoiding W island formation. 
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• Increased flow rates to minimize resident time and product gases. 

Furthermore, low growth rate implies longer processing time,- which again causes 
loss of selectivity. Therefore although selective tungsten deposition is desirable because 
it avoids the additional step of etch back, it is still not the choice for manufacturability. 

From the discussion above, it is clear that still blanket W deposition processes has 
to be used for depositing W in trenches of ICs. This allows much more relaxed control 
and also the throughput is large. However, the blanket tungsten deposition suffers from a 
problem of void formation[l 7], which may affect the reliability of interconnects in ICs. 

The experimental work of Chen et al.[18] illustrates the deposition of blanket 
tungsten layer followed by etch back. Depositions are performed on both patterned and 
vmpattemed wafers to study the plug formation and to characterize the film properties. 
Excellent step coverage is reported for gas ratio 1:2 of SfflU to WFe. The etch back is 
done by SFe and N 2 . 

Apart from experimental studies, literature provides many models for filling a 
trench during blanket deposition. Jain et al.[9] provide a comparison of a diffusion 
reaction model and a ballistic transport reaction model and find a particle based line of 
sight model to give better predictions. However, this approach is limited to regimes 
where assumption of Knudsen flow is valid. 

Hasper et al.[19] presented a model to calculate the step coverage of blanket 
tungsten low-pressure chemical vapour deposition (LPCVD) by hydrogen reduction as 
well as silicon reduction. For hydrogen reduction, a good agreement between 
experimental and model step coverage is obtained. HF is found to have no influence on 
step coverage. The step coverage, in high aspect ratio (2|im x lOjam) trenches is 
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enhanced by high WFe concentration and low growth rate. They obtained a step 
coverage of 95% in trenches with aspect ratio 5 uniformly over a 6 inch wafer. Schmitz 
and Hasper[17] also studied the step coverage of blanket W CVD for aspect ratios 0.5, 
1.0 and 2.0. The flow rate of WFg is 20 seem., and that of H 2 is 500 seem., at a 
temperature of 430 °C. They found that the size of void increases with increase in aspect 
ratio of the trench. The results are also experimentally verified. Egashira et al.[20] also 
presented a model on blanket CVD of tungsten. Here they have described the reaction 
mechanism of WSix blanket CVD by WF 6 /Si 2 H 6 . They have also predicted the reactor 
performance based on the data of a tubular reactor. The simulation is aimed at the 
optimization of a newly developed reactor. 

Another model by Klejin et al.[21] studies the influence of WFe concentration on 
growth rate in W LPCVD by hydrogen reduction, both experimentally and 
mathematically. Model predictions are in good agreement with experimental growth 
rates. Above a certain critical partial pressure of WFg, growth rate is found to be 
independent of partial pressure of WFe. Thermal diffusion in cold wall reactor is found to 
be very important. 

Most models reported in literature are two-dimensional. A two dimensional 
model is sufficient to provide insight into reasons for void formation. However, a three- 
dimensional model is presented by Bar and Lorenz[22], which is more realistic and it can 
be applicable to arbitrary shapes. 

Yim and Park[4] studied the numerical solution for high aspect ratio trench 
models to study the step coverage in LPCVD process. They use Monte Carlo simulation 
and string algorithm. Sticking coefficient, surface diffusion coefficient, direct deposition 
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and various re-emission mechamsms are included in this model. It is found that the step 
coverage becomes poor as the temperature and aspect ratio increases. Experimental 
results do agree with simulated ones. Finally, they concluded that low temperature 
operation below 280 °C is proposed as a way to deposit the conformal tungsten film with 
good step coverage. 

In another Monte Carlo simulation model for W LPCVD on silicon substrate by 
hydrogen reduction, Kim et al.[23] showed the non-lineffl: dependence of WFe 
concentration on growth rate. The study adopts a simple-surface reaction model 
assuming adsorption of WF6 only on trench walls. They showed that growth rates 
increase linearly. 

From this review, it is clear that blanket W deposition is more manufacturable 
than selective W deposition. Machines are now available, in which both deposition and 
subsequent etch back process can be integrated. But, both experimental and model 
studies demonstrate the formation of a void during the filling process of the trench. The 
void formation becomes more pronounced as the deposition temperature and aspect ratio 
of the trench increases, which is eventually a reliability issue. 

1 .3 OBJECTIVES OF THE PRESENT STUDY 

A silicon wafer contains a number of devices, electrically connected by 
interconnections. These multiple layers of interconnects are connected by tungsten plugs. 
The shapes of these plugs are, for example, rectangular trenches having width less than a 
pm and length of a few pm. For making contacts of tungsten in these trenches, WF(, is 
passed over the wafer along with hydrogen. These gases diffuse into the trench, where 
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they react at the walls of the trench. As a result tungsten is deposited on the walls and HP 
diffuses out. Because these source gases have to diffuse throughout the depth of the 
trench, concentration gradient may develop, with higher concentration at the feature 
mouth. Since the concentration of the reacting source gases is more at the mouth, the 
reaction rate is also faster there than at the bottom of the trench, that is, the rate of growth 
of tungsten layer is more at the mouth. This leads to the closure of the mouth before 
complete filling of the trench, resulting in formation of a void in the trench. This feature 
can adversely affect device reliability and device performance. Schematically, this, 
situation is shown in Fig.l. 

The objective of the present work is to ensure voidless deposition of tungsten in 
the rectangular trenches of the interconnects of ICs. Optimizing process parameters that 
give the best step coverage can do this. In other words, we have to proceed with 
deposition in kinetic control regime, where concentration gradients are minimized. This 
is the reason why we get better step coverage at low temperature. But at low temperature 
the deposition time is very large, which affects the throughput. 



Fig. 1 : Schematic of a void formation in blanket W CVD. 
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Obviously, it will be much better if W CVD can be carried out at higher 
temperature and (atmospheric) pressure without void formation. Thus the goal of this 
work is to explore new ways in which this can be accomplished. One way would be to 
impose a temperature gradient along the vertical direction of the trench, so that, its 
bottom will be hotter than the mouth. As a result, W will preferentially deposit at the 
base of the mouth and the trench would fill inside-out. 


12 



CHAPTER 2 


CHEMICAL VAPOUR DEPOSTTTON : A TWO 
DIMENSIONAL MODET. 

2.1 THE MODEL 

As discussed in the preceding chapter the process of filling trenches of tungsten 
in interconnects of ICs is performed by chemical vapour deposition (CVD) process. 
This chapter deals with two-dimensional modeling of a trench filling process by CVD. 

In VLSI device fabrication process, each device contains submicron trenches/vias. 
The shape of these features is taken as a rectangular section. A cross section of the 
trench, infinitely long in the third dimension is shown in Fig. 1 . 

^ dielectric 


Si or metal 



Fig.l : Schematic of a cross-section of a trench 
We consider only one half of the trench as the computational domain for mass 
transport calculations, allowing for symmetry about x=0 (see Fig.2). The initial length 
and half-width of the trench is taken as Lq and Wq, respectively. 
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Fig.2: The computational domain 

Deposition of tungsten by CVD using WFe as the source gas proceeds according 
to the following heterogeneous chemical reaction: 

WF6(g) + SHaCg) = W(s) + 6HF(g) 

The product gas HF diffuses out of the trench. The rate of tungsten deposition by 
H 2 reduction of WFe has been reported as[9]: 

«=2.2649 xlO-'exp po^x^ 

where kf= 0.75006 (dyne/cm^)’’, P is the partial pressure of WF6 at the mouth 

of the via. For any significant value of P ^p^ , the term in the square bracket in equation 
(1) goes to unity. 

An expression for molar flux is determined from ‘dusty’ gas model[27]. This 
model takes into account ordinary diffusion, viscous flow and Knudsen flow and has 
been successfully applied before[28]. The derivation of flux expression starts from 
kinetic theory of gases. Here the multicomponent diffusion of gases occurs by binary 
collision between gas molecules. The solid in this model is also considered a heavy, 
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immobile gas. Knudsen diffusion is included assuming binary collision between gas 
molecules and the dust gas. Furthermore, viscous flow is also considered. 

Neglecting thermal diffusion, the flux expression can be written as follows: 


i5:^[x,N,-X,Nj.i = -VX,-X. 


1 W- 
P SpDi, 


VP 


( 2 ) 


for ij =1,2,3 

Subscripts 1,2,3 indicate H 2 , WFe and HF respectively. The quantities Dik, Oy and 
T] are estimated by the formulae given in Appendix I. In writing this equation, we have 
used the fact that for a flow through a region between two parallel plates, the average 
velocity in direction of pressure gradient is: 


V =■ 

avg. 


3^ 


VP 


( 3 ) 


Under assumption of quasi-steady state, the mass balance for H 2 is 


V.N, = 0 (4) 

Also, from the stochiometry of the hmgsten deposition reaction, 

N2=-N, (5) 

' 3 ' 

N3=-2Ni (6) 

In addition, we have ideal gas law 

P = cRgT (2) 

And finally, 




( 8 ) 
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Expanding eq. (2), we can write three equations in terms of molar fluxes of H2, 
making use of eq. (5) to (8) at the same time. 

For i = 1, 




•cVX, - 

R.T 


or, N, =a, VX, +p, VP ■ (10) 

For i = 2, 



or, N, = a, VX, + + P3 VP (12) 

where, all a’s and P’s are listed in Appendix II. 

Now, when eq. (10) to (12) are substituted in eq. (4), we get three partial 
differential equations containing variables Xi, X2 and P. These are solved using the 
following boundary conditions. 

If L(t) is the depth of the trench any time, t, and the deposition profile is as given 
in Fig. 3, then 
at, x=0, N,.n, =0. 

at, y=L(t), that means, at the trench entrance 
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Xi=X°, 


X2=X^ 

p=p°. 

and at the tungsten moving boundary, N, =3iR 



5 ^ 

Fig.3: Schematic showing boundary conditions 

2.2 METHOD OF DISCRETIZATION AND SOLUTION 
PROCEDURE 

The computational domain at t=0 is shown in Fig. 2. The partial differential eq. 
(2) is solved by a standard finite difference method. This computational domain is now 
discretized into a number of mesh points. Under the assumption of quasi-steady state, the 
finite difference form of eq. (4) at any internal point, shown in Fig. 4, in the mesh can be 
written as follows. 

[n: -N;:] (y, - y„)+ [n; - x.)= 0 (15) 
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Where, N, - N^x + N^y and superscripts represent the location in reference, to 

Fig. 4. Point N, E, S and W are points on north, east, south and west, respectively, of 
point P and n, e, s and w are respective mid points. 


N 



Fig.4: Discretization at a node point 

There are three expressions for the flux, eq. (10)-(12), which Mien substituted in 
eq. (15) giving three finite difference equations in Xi, X 2 and P at any one node point. 

As an example, when eq. (10) is substituted in eq. (15), an expression for N® becomes 



where xe and xp are x-locations of points E and P, respectively, aie is calculated as: 

«j_ 

^ , in which aiE and aip are values of ai at points E and P, respectively. 

In this manner, by substituting eq. (10)-(12) in eq. (15), three discretized 
equations are obtained. At any particular grid point of the mesh, the three discretized 
equations contain variables Xi, X 2 and P. The complete set of discretized equations are 
given in Appendix-Ill, which also includes boundary node points. 
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Clearly, the quantities a’s and P’s used to discretize the PDE’s are dependent on 
the variables Xi, X 2 and P themselves. So in order to compute a’s and p’s we use a set 
of guessed values for Xi, X 2 and P; then Xj, X 2 and P are computed at each node point. 
If the guessed values and the computed values match within 10'*, then convergence is 
achieved. Otherwise, the computed values for Xi, X 2 and P are now treated as guessed 
values and the procedure is repeated. 

These finite difference discretized equations were solved by using the subroutines 
BANDEC and BANBKS[30]. These subroutines employ LU decomposition method for 
solving the linearised equations. 


2.3 GROWTH AND BOUNDARY MOVEMENT 

The growth rate, that is, rate of change in tungsten film thickness at any point on 
the wall of the via is written as: 


Mw PX, f l + k,P°xf 

dt ' P“X“ l + kfPX^ J 


(17) 


Initially, W = Wo at t = 0. 

Once the concentrations (Xi, X 2 and P) are computed at all nodes, the growth at 
any node point is determined according to the reaction rate in eq. (1). Thus the boundary 
in a direction normal to itself moves by an amount 


p„ P”x; I I + k,PX, 


At this stage, the growth algorithm is best discussed by considering growth at t=0 


during a period At; in all subsequent time steps, the algorithm remains the same. 
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Consider the mesh at t=0 imposed on left half of the trench, the solid lines in Fig.5, where 
at each node point we have already computed the concentrations. In the mesh, we have 
to determine the movement of the boundary points labeled “a”, “b” and “c”. The points 
“c” only move up, in -Y direction and points “b” move only in -X direction. Whereas, 
point “a”, the left most point on the trench mouth, where boundary condition specifies 

concentrations X ° , X ” and P ° , moves in both -X (to a') and -Y (to a") directions. 



c 

Fig.5 : Point a, b and c on the mesh, equally spaced in 
X-direction and logarithmically in Y-direction 
Since the concentrations are known at “a” a priori, the time step At is computed as 
time taken to move to “a'”. Thus in time At, the boundary at point “a” moves to “a'”. In 
the same duration, by an equal amount, W grows on the surface of the wafer. This fact is 
represented by increasing the depth of trench and adding an additional row of mesh 
points, seen as dotted lines in Fig. 5. These added mesh points represent the trench 
entrance for concentration calculations in next time step, after incrementing time by At. 
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Thus the point “a” also moves to point “a"”, which in immediately following iteration 
plays the same role as “a” in this iteration. 

Notice that the mesh in X direction is equally spaced. However, in the other 
direction, near the trench entrance, the mesh spacing increases logarithmically and then 
becomes equally spaced. If the mesh in Y direction were equally spaced, because 
Lo»Wo, the mesh spacing in Y direction will be much larger than that in X direction. 
And because in Y direction, point “a” grows to “a"” by an amount equal to spacing in X 
direction (that is distance between “a” and “a'”), the vertical spacing for the newly added 
row will be much less than the adjacent spacing between points “a” and neighboring . 
point “b”. In order to represent the jSnite difference approximation accurately, we require 
that the mesh spacing in a direction should not increase by more than a factor of 1.2. 
Thus, the vertical direction spacing near the trench entrance is increased logarithmically, 
keeping spacing between “a” and adjacent “b” within a factor 1.2 of X direction spacing. 

So far we have discussed the computation of time step and movement of point, 
“a”. Now we establish general strategy for movement of other boundary points labeled 
“b” and “c”. Each of these points moves according to eq. (18). In time At, there is no 
guarantee that a point will move to another node point. For example, in Fig.6, a point “b” 
moves the boundary by an amount Ax, to a node point in between “b”’ and “b"”. 

Ax ^ 

• • • 

b b' b" 

Fig.6: Depiction of boundary movement along X-direction. 
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However, we require that for the purpose of calculation of concentrations in the 
next time step, the boundary lies on a node point. This difSculty is avoided by keeping 
track of the actual movement of boundary at any node point, but for the purpose of 
calculation of concentrations, placing the boimdary point on the nearest node point. Thus 
in the example above, while the actual displacement Ax is recorded and used in all 
subsequent calculation of boundary movement, the concentration calculations are 
performed by placing the boundary on the nearer of the points “b'” and “b"” to determine 
Ax. 

Further, growth of points labeled “c” is computed. We find that the concentration 
gradients in X direction are not sufficient to cause points “c” to grow by an unequal 
amount. Thus to simplify the problem, the growth of the lower boundary (points “c”) is 
regarded as uniform. After the growth of the boundary to new points, these new 
boundary points are labeled as points “c” in the next iteration. 

Once the growth of the lower boundary is computed, the nodes lying fully in W 
(that is other than the new boundary nodes) are rendered useless in future time steps. 
Also the boundary points “b” consumed due to movement of lower boundary are also 
removed from the reckoning; then movement of boundary at remaining “b” points is 
computed. According to the amount of growth, the new node points obtained as 
boundary "are labeled as points “b” for the next iteration. We ensure that this new 
boundary at adjacent points “b” does not differ by one node point in X-direction; 
otherwise the finite difference equation in Appendix III pertinent to this boundary will 
not apply. In the event, this requirement is violated we start again with a finer mesh. 
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List of Symbols used in Chapter 2 

local gas concentration in moles per unit volume 

binary diffusion coefficient between species i and j 

ICnudsen diffusion coefficient of species i 

first order heterogenous rate constant 

Boltzman constant 

constant used in rate expression 

initial depth of the trench 

depth of trench at any time t 

molecular weight of species i 

molecular weight of tungsten 

molar flux of species i 

total pressure 

partial pressure of species i 

partial pressure of species i at the trench mouth 

rate of heterogeneous chemical reaction 

universal gas constant 

time 

temperature 

average velocity in direction of pressure gradient 

local half-width of the trench at any y 
initial half-width of the trench 
mole fraction of species i 
viscosity of gas mixture 
density of tungsten hexafluoride 
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APPENDIX I 


The knudsen diffiision coefficient is given by the expression 


2 8RJ 

3 \ m,-7r 


where nij is the molecular weight gas species i. 

From kinetic theory of gases the viscosity of gas species i is given as 


Tii =2.6693x10 


-5 




of 

where rji is in gm cm‘’s'', mi is in gm, T in K and Lennard-Jones parameter cti is 
0 22 

in A . Q is the collision integral tabulated as a function of dimensionless temperature 
T* = ksT/ei,, where ei is another Lennard-Jones parameter. 

The viscosity of gas mixture is given by Wilke’s semi-empirical formula: 


XjTlj 

j = l 




where O.-j 

" 4 


1 + 






v®iy 


iH-i^ 


Again firom kinetic theory of gases an expression for binary diffusion coefficient 


is given by; 


a. = 1.8583x10 


-3 



PcyfjQ‘'T.. 
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/ 

where Sjj is in cm^/s, T in K, P in atm,CTjj = 

V 


wheree.. = ^e,.G^. . 



APPENDIX TT 


The expressions for a ’s and p ’s are as follows: 
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APPENDIX TTT 


The expression for rate of tungsten deposition from WFe by hydrogen reduction is 
given as a function of Xi, X 2 and P. 


9?(Xi,X2,P) = kP°-^X-°-^ 


P°X,° 


l + k^P^X, 


Ov o^ 


l + kfPXj j 

Using multi-variable Taylor’s expansion series, all the three variables can be 
separated. 


iR(X,,X2,P) = iR(X,®, X2^ P®) + (X, - X,8) 


5R 


aXj 


+ (X2-X2®) 


SR 


sx. 


+ (P-P®) 


SP 


Where Xi® , X 2 ® and P® are the guess values at every node point. 

For convenience, we can express three times of rate as: 

39?(Xi,X2,P) = Co + ciXi + C 2 X 2 + C 3 P 

Now using this expression for iR(Xi,X 2 ,P) we can write the discretized equations 
for every node of the mesh. 


For a node present on ‘b’ of Fig.5: 

Considering the fluxes only in X-direction, we can write 

[N: -91(X,,X2,P)](ys-yn) = 0 

Using the flux expression given in equation (16), the above equation can be rewritten for 
all the three components. 
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In a similar fashion, the discretized equations in other regions of the trench can be 
obtained. 
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For a node present on ‘c’ of Fi P.'S r 
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For internal nodes present in Fig.5! 
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For the nodes present on right boimdary of Fip 
[N: -N:](ys-yn) + [N; -N;](xw-Xe) = 0 
Applying image point technique, N* = -N* 
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CHAPTER ^ 


■RESULTS AND DISCIISSTON 

The computed growth profiles for blanket tungsten CVD are discussed in this 
chapter. The values of parameters used in calculating the growth profile of W deposition 
process is given in Table 1 . While Lo is always 2pm, we change the aspect ratio, which 
is the depth divided by the half-width of the trench, by varying Wo between 0.2pm and 
0.5pm. 


Table 1: Values of parameters used in computing results 


Density of tungsten, pw, gm/cc 

19.22 

Molecular weight of tungsten, Mw, g/mole 

183.85 

Constant used in rate expression for tungsten deposition, kf, (dynes/cm"^)'* 

0.75006 

Mole fraction of H 2 , X ^ 

0.9615 

Mole fraction of WFa, X 2 

0.0384 

Mole fraction of HF, X 3 

0.0 

Total pressure, P'^, dynes/cm^ 

1333.2 

Depth of the trench, Lo, pm 

2.0 


For calculating binary diffusion coefficient and viscosity (as given in Appendix I), 

11 22 

we need Lennard-Jones parameters and values of collision integrals, Q and Q . The 
lennard-Jones parameters are provided in Table 2[11], whereas the values of Q and Q 
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are interpolated at the corresponding T (-keT/e) by using standard table[30] for collision 
integral. 


Table 2: Lennard-Jones parameters 


0 

CT, H 2 , A 

2.827 

e/ke, H 2 , K 

59.7 

0 . 

a, WFfi, A 

5.21 

e/ks, WFe, K 

338 

0 

a, HF, A 

3.15 

E/ks, HF, K 

330 


First sets of results are on comparison of mesh. For this we have taken a trench of 
width Wo=0.35pm at 773 K. the computational domain is divided into 100 interval both 
in X and Y direction. The progression of deposition of time is depicted in Fig.l. In 
plotting this and subsequent similar figures, we have shifted the origin to center of the 
trench base (see Fig.2). 

Also, the various deposition profiles are for increasing time, x, expressed as actual 
time divided by total processing time; the process is deemed oyer when W/W o is zero at 
L/Lo=l . The plotted profiles are for x ranging from 0 to 1 , in intervals of 0. 1 . 

The results in Fig.l demonstrate that in the initial deposition period, the trench 
fills uniformly. However, as the trench width becomes small, process becomes mass 
transport limited, and a concentration gradient sets in. Hence the trench entrance fills 
faster, eventually leading to a void. 
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W/Wq 


Fig. 1 : Isothermal deposition profile at 773 K for Wg -0.35pm and 
mesh size 101x101; total processing time is 40 seconds 
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The formation of a void is consistent with both experimental observation[ 17] and 
modeling of it[9,23]. In order to determine the mesh sensitivity in onr calculations, we 
next divided the domain in 50 intervals, instead of 100, in both X and Y directions. 



WAVo 


Fig.2: Schematic of a trench corresponding to Fig.l . 


The formation of a void is consistent with both experimental observation[17] and 
modeling of it[9,23]. In order to determine the mesh sensitivity in our calculations, we 
next divided the domain in 50 intervals, instead of 100, in both X and Y directions. The 
result for this case, in Fig.3, is compared with that in Fig.l. We find no significant 
difference between the two; the normalized void volume, defined as volume of the void 
divided by initial trench volume, is 7.0% and 6.7% for cases in Fig.l and Fig.3, 
respectively. Since this difference is not significant, and that it took 4 hrs and 30 minutes 
of real time on a pentium II machine with 50x50 mesh, while 3 days and 18 hrs in &e 
case of 100x100 mesh, in all subsequent calculations, the mesh used is the same as in 
Fig.3. 
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Fig.3: Isothermal deposition profile at 773 K for Wg -0.35p,m 
and for a mesh 51x51; total processing time is 40 seconds 





3.1 DEPOSITION AT 773 K 

The first few computations are for W deposition under isothermal conditions. In 
addition to results presented in Fig.3, at a deposition temperature of 773 K, which is 
regarded as an upper limit of deposition temperature, the simulation is performed for 
three additional trench widths 0.2pm, 0.25pm and 0.5pm. 

Fig.4 depicts the growth profile for a trench 0.5pm wide, where a void is found to 
be formed having a normalized void volume 3.4%. Thus the normalized void volume 
increases from 3.4 to 6.7 as trench width decreases from 0.5pm to 0.35pm (see Fig.3). A 
further increase in the void formation is noticed when Wo decreases to 0.25pm and 
0.2pm. Fig.5 and 6 show the simulated growth profiles for trench width 0.25pm and 
0.2pm. The corresponding normalized void volumes are 12.2% and 17.3%, respectively. 

Thus, at a deposition temperature 773 K, a significant void forms in all the cases. 
This is due to the concentration gradient that develops along the vertical direction of the 
trench. Since, the reaction rate is dependent on concentration of reacting species, the 
mouth of the trench fills faster. Thus, once the trench entrance is closed, a void forms. 

Also, the normalized void volume increases, as the aspect ratio of the trench 
increases. Clearly, as the trench width decreases, more resistance is offered to the 
diffusion of reacting species, which, in turn, shifts the deposition process to mass transfer 
control. As a result, stiffer concentration gradient is set in vertical direction. This leads 
to a more rapid mouth closure in high aspect ratio trenches. 

In prelude to this work, we have suggested that imposing a temperature gradient 
may alleviate the problem of void formation. If due to concentration gradient, the 
reaction rate is faster at the trench entrance than at the bottom of the trench, we impose a 
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Fig.4: Isothermal deposition profile at 773 K for =0.5|a.m; 
total processing time is 57 seconds 
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temperature gradient along the depth of the trench with higher temperature at the bottom 
than at the entrance. This would compensate for lower deposition rate at the bottom of the 
trench. Here, in all four cases, we kept the entrance of the trench at 773 K, but linearly 
increase the temperature to the bottom by 3 K; no temperature gradient is imposed in X 
direction. After every successive growth step, the depth of the trench increases and again 
the mouth temperature is fixed at 773 K. The simulated profiles are shown in Fig.7 
through 10. These figures depict that, in no case, by imposing a temperature gradient 
avoids void formation. However, we are able to decrease the normalized trench volume 
in all the cases. The normalized trench volumes are decreased from 3.5, 6.7, 12.2 and 
17.3% to corresponding values 2.1, 5.3, 10.7 and 16.1% for trench 0.5imi, 0.35iim, 
0.25 pm and 0.2pm wide, respectively. Thus from the discussions so far, it is obvious 
that 773 K is not an ideal deposition temperature for W CVD. It would be possible, to 
impose a steeper temperature gradient in order to completely fill the void. But such 
gradients may be difficult to obtain. Thus we fix the maximum difference of temperature 
(AT) between entrance and bottom of trench to 3 K, and examine void formation at a 
lower temperature. 

3.2 DEPOSITION AT 673 K 

Few simulations are performed at a deposition temperature 673 K. In Fig.l 1 we 
represent the isothermal profile for a trench of width 0.5pm. Now the normalized void 
volume is 0.5%, much lower than 3.4% at 773 K. Similarly, in Fig.12-14, we show the 
growth profile at 673 K for trench widths 0.35pm, 0.25pm and 0.2pm, respectively. The 
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W/Wo 


Fig. 8: Modified deposition profile of Fig.3 under non-isothermal 
condition; total processing time is 40 seconds 
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Fig.9: Modified deposition profile of Fig.5 under non-isothermal 
condition; total processing time is 28 seconds 
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Fig. 10: Modified deposition profile of Fig.6 under non-isoliiermal 
condition; total processing time is 23 seconds 
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Fig. 1 1 : Isothermal deposition profile at 673 K for -0.5pm; 
total deposition time is 308 seconds 
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Fig. 12: Isothermal deposition profile at 673 K for -0.35pm; 
total deposition time is 215 seconds 
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Fig. 13: Isothermal deposition profile at 673 K for =0.25pm; 
total deposition time is 154 seconds 
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Fig. 14: Isothennal deposition profile at 673 K for -0.2pm; 
total deposition time is 123 seconds 
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corresponding normalized void volumes are 1.2%, 2.4% and 3.5%. These values are also 
much lower than the void volume at 773 K. 

However, also at 673 K, a void forms in all cases discussed here. However, the 
normalized void volume is comparatively lesser, when compared to the corresponding 
cases at 773 JC. That means, for a particular aspect ratio, the normalized trench volume 
increases with an increase in temperature, dearly, with increase in deposition 
temperature, the reaction rate increases exponentially while the diffusion constants 
increase by a smaller power of temperature. Thus the reaction kinetics becomes faster 
allowing the deposition process to shift further to mass transport controlled regime; as a 
result, keyholes are formed which are larger in volume than those formed at lower 
deposition temperatures. 

Here at 673 K also, a controlled thermal gradient is applied in the vertical 
direction of the trench (as in the cases for 773 K). As depicted in Fig. 15, for a trench 
width of O.Spm, only 1 K temperature difference between entrance and the bottom of the 
trench is sufficient to ensure complete filling of the trench. But, for the trench width 
0.35p,m (see Fig.l6), 2 K temperature difference is needed for void-ftee deposition. 
However, for trench widths 0.25pm and 0.2pm, as shown in Fig.l7 and Fig.18, 
respectively, we need a temperature difference of 3 K for complete filling. 

The results clearly demonstrate that the imposition of controlled thermal gradient 
can be employed to completely fill the trenches. 
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0.2 0.4 0.6 0.8 1.0 

W/Wo 


Fig. 1 5: Modified deposition profile of Fig. 1 1 after application of AT-1 K; 
total processing time is 307 seconds 
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Fig. 17: Modified deposition profile of Fig.l3 after application of AT 3 K; 

total processing time is ,154 seconds 
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Fig. 18: Modified deposition profile of Fig.l4 after application of AT-3 K; 
total processing time is 123 seconds 





3.3 DEPOSITION AT 573 K 


Indeed, just as going from 773 K to 673 K lead to reduction in normalized void 
volume, a further reduction in temperature for isothermal CVD could also ensure 
complete filling of the trench. This is demonstrated by reducing the deposition 
temperature to 573 K. 

In all four cases corresponding to Wo=0.5, 0.35, 0.25 and 0.2pm, a completely 
filled trench is achieved without any application of thermal gradient (see Fig. 19-22). 
Thus it has been necessary to conduct CVD at a lower temperature to bring deposition 
under kinetic control and hence uniform deposition. However, this requires longer 
processing time, and hence longer through put. 

Thus we have shown for the first time that by apphcation of temperature gradient, 
it will be possible to fill a trench at a higher temperature. This would result in an 
increased through put. 

3.3 THE THERMAL GRADIENT 

In this work so far, we have proposed a maximum gradient of 1 .5K/pm. We now 
discuss the feasibility of supporting such a temperature gradient. 

First we assume that the gases in the trench are in thermal equilibrium with the 
neighboring dielectric. If the wafer loses heat by convection from the top surface, then 
the heat transfer at t=0 is defined according to, 

Kdid^ = l»(Ts-298). 

^0 
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Fig. 19: Isothermal deposition profile at 573 K for -0.5pm; 
total deposition time is 3004 seconds 
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Wo=0.2nm 
T=573 K 



Fig. 22: Isothermal deposition profile at 573 K for WQ=0.2pm; 
total deposition time is 1207 seconds 
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where Kdiei is the thermal conductivity of the dielectric, AT is the temperature difference, 
h is heat transfer coefficient and Tj is the temperature of the wafer surfece. For example, 
if we take the heat transfer coefficient as 500W/m^K (upper limit for forced convection), 
and surface Ts-298 as 500K, we find to support IK temperature difference over Lo, die 
dielectric should have a thermal conductivity of 0.5 W/mK. 

The most common insulation/dielectric used in silicon technology is Si02. The 
thermal conductivity of this material in various forms ranges between 0.6-1. 5 W/mK[35]. 
Thus, the conventional ceramic is going to be just at the upper limit of required value of 
thermal conductivity. 

As the aspect ratio of the trenches becomes more severe, polyimide film.s become 
attractive for interlayer dielectric. The polyimide materials have conductivity in range 
0.2-0.3 W/mK[36]. 

Currently, we are working with 0.18pm design rules, metal pitch 0.5pm and 
aspect ratios of trenches as 2.5. In such cases, conventional dielectric are being replaced 
by much lower dielectric constant materials such as Fluoropolymer films[37]. Examples 
in this class includes materials such as Parylene N &F, teflon, poly(napthalene) PNT-N 
and poly(fluorinated napthalene) PNT-F. Most measurements on these materials are on 
dielectric constant, which is lower than that of Si 02 ; less is known about their thermal 
conductivity. We will have to explore the thermal conductivity values for these 
materials, and then the application of required thermal gradient may become more 
convenient. 
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CHAPTF.R 4 


CONCLTJSTON 


In this work, by a two-dimensional model, we have studied in detail blanket 
tungsten deposition by hydrogen reduction of tungsten hexafluoride in a long rectangular 
trench. The effect of imposing a controlled thermal gradient, along the vertical direction 
of the trench is also investigated. The temperature gradient is such that, higher 
temperature is maintained at the base of the trench, compared to the entrance of the 
trench. Then these temperatures are input to the mass transport and chemical reaction 
model. 

A sophisticated mathematical model was presented here, which deals with the 
diffusion and chemical reaction in a CVD reactor. The model assumes ordinary 
diffusion, viscous flow and Knudsen flow in a multicomponent system. 

The dusty gas based flux expression[27] is used to generate a system of partial 
differential equations, containing variables: mole fraction of hydrogen (Xi), mole fraction 
of tungsten hexafluoride (X 2 ) and total pressure (P). The set of partial differential 
equations (PDE’s) with appropriate boundary conditions is solved, using a control- 
volume-based finite difference method. The equations are discretized in a non-uniform 
grid and are solved by LU decomposition method[29]. The grid independence is verified 
by using 51x51 and 101x101 grids. Very little difference was found between the two 
results. 

From the results of deposition of tungsten under isothermal condition, we find 
that at a lower deposition temperature (573 K), a trench 2pm deep and of width 0.5pm, 
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temperatures (673 and 773 K), incomplete filling of the trench is shown, leading to the, 
formation of a void inside the trench. At a particular deposition temperature, the size of 
the void increases as the trench aspect ratio increases. Moreover, for a particular aspect 
ratio, the void volume becomes larger with increasing temperature. 

Thus it is concluded that a lower deposition temperature and a larger initial trench 
width promotes conformal tungsten deposition with greater step coverage. However, this 
prolongs the deposition time to fill a trench, which affects the throughput. Obviously, it 
will be much better if void free W CVD can be carried out at higher deposition 
temperatures. 

Growth of tungsten without a void formation can be achieved by imposing a 
thermal gradient. This idea has never been tried before. 

In our two-dimensional model, we show that the W deposition, imder an imposed 
linear temperature gradient results in a smaller normalized void volume. For deposition 
at 673 K, we show that by applying a gradient of less than 1.5 K/pm, a trench can be 
completely filled, while corresponding isothermal processing lead to a void. Thus the 
model presented here could help the equipment designers and process engineers to 
establish a reliable process for blanket tungsten deposition. 

Finally, in order to attain a thermal gradient of 1.5 K/pm, the conventional SiOa 
dielectric will not be adequate. However, replacing Si02 with polyamide and fluoro- 
polyamide could make such thermal gradient feasible. 
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